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ABSTRACT 

The  mathematical  model  used  in  previous  Davidson  Laboratory  adapta¬ 
tions  of  linearized  unsteady  lifting  surface  theory  to  marine  propellers 
has  been  revised  by  removing  the  so-called  "staircase"  approximation  of 
the  blade  wake  and  replacing  it  by  an  "exact"  heiicoidal  blade  w,  ke.  A 
new  numerical  procedure  and  program  based  on  the  present  model  has  been 
developed  to  evaluate  the  steady  and  unsteady  blade  loading  distributions, 
which  are  used  to  determine  the  bearing  forces  and  moments.  Systematic 
calculations  of  these  forces  and  moments  for  a  series  of  propellers  show 
better  agreement  on  the  whole  with  experimental  measurements  than  did 
the  earlier  calculations  for  the  same  series.  In  addition,  the  chordwise 
loading  distributions  are  much  smoother  than  obtair.eJ  previously.  However, 
the  quantitative  improvement  must  be  weighed  against  the  considerable 
increase  In  computer  time  over  the  old  method. 


Keywords 

Hydrodynamics 

Unsteady  Theory  for  Marine  Propellers 


Mi 


R-1509 


TABLE  OF  CONTENTS 

ABSTRACT  . .  Ml 

NOMENCLATURE .  vll 

INTRODUCTION . . .  I 

LINEARIZED  UNSTEADY  LIFTING -SURFACE  THEORY  .  3 

DEVELOPMENT  OF  THE  KERNEL  FUNCTION  .  4 

REDUCTION  TO  A  ONE-DIMENSIONAL  INTEGRAL  EQUATION  .  10 

PROPELLER  LOADING  AND  RESULTING  FORCES  AND  MOMENTS  .  IS 

Loading .  16 

Propeller-Generated  Forces  and  Moments  .  16 

Blade  Bending  Moments .  17 

NUMERICAL  RESULTS  .  18 

CONCLUSIONS  .  23 

REFERENCES .  25 

FIGURES 

APPENDIX  A:  Singularity  of  the  Kernel  Function  at  p  ■  r 

APPENDIX  B:  Evaluation  of  the  0  -  and  $  -integrals 

a 

APPENDIX  C:  Evaluation  of  the  Integrand  of  the  x-lntegral  at 
the  Singularity  \  ■  0 

APPENDIX  D:  Evaluation  by  the  Lagrange  Interpolation  Method  of 
the  P-lntegratlon  In  the  Region  of  the  Singularity 


v 


R-1509 


a 


b(r) 


Fx,y,z 

f(r,  s) 


g(x) 

•m<> 

,(m)(  } 

.po 

<,K>m 

I 

V) 

K 

*m<> 

fc 

k(«".n) 

L* 

L(r) 

j e 

Hb 


NOBfcttSLATURE 

Q/U 

function  defined  in  Eq.  (28) 

expanded  semichord  length  at  each  radius,  ft 

propeller-induced  forces  in  x,y,z  direction 

camberline  ordinates  from  face  pitch  line  at  each 
radial  position 

function  defined  in  Eq.  (28) 

modified  Bessel  function  cf  first  kind,  or  order  m 

defined  in  Appendix  B 

defined  in  Appendix  B 

defined  in  Eq.  (28) 

index 

Bessel  function  of  order  m 

kernel  of  integral  equation 

kernel  after  chordwise  integrations 

modified  Bessel  function  of  second  kind,  of  order  m 

variable  of  integration 

influence  functions  (Eq.  40) 

loading  distribution  in  lb/ft^ 

spanwise  loading  distribution  in  lb/ft 

spanwise  loading  coefficients  of  Bimbaum  distribution 

integer  multiple 

blade  bending  moment  about  face  pitch  line 
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final  value  of  m  In  summation 
Index  of  summation 
order  of  lift  operator 
number  of  blades 
blade  Index 
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pressure 
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Fourier  coefficients  of  velocity  distribution 
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velocity  due  to  flow- incidence  angles 
velocity  distribution  normal  to  propeller 
longitudinal  ordinate  of  control  point 
cylindrical  coordinate  system  of  control  point 
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e(ii) 


ep(r) 


A<S)(x) 
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transverse  cartesian  ordinate 

vertical  cartesian  ordinate 

angle  in  Eg.  (18) 

half-length  of  radial  x  trip 

small  length 

cho/dwisu  modes  (Eq.  21) 

angu >r  ordinate  of  loading  point 

ang  * posit 'xt  rf  loading  point  *ith  respect  to 
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proje  *.ed  blade  **•?*»  &rd  lengtn,  -ad  Ians 
(n-i) ,  n  *  i , 

geometric  ;**ch  angle  at  each  radial  position 

angular  ohordwls'?  location  of  loading  point 

defined  in  Appendix  B 

defined  in  Appendix  B 

order  of  harmonic  of  loading  distribution 

variable  of  integration 

longitudinal  ordinate  of  loading  point 

cylindrical  coordinate  system  of  loading  point 

radial  ordinate  of  loading  point 

mass  density  of  fluid 

angular  measure  of  skewness 

or 

Cesaro  sum  of  curtailed  Birnbaum  series 
variable  of  integration 
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INTRODUCTION 

In  the  course  of  a  series  of  Investigations  concerned  with  adapting 
linearized  unsteady  lifting  surface  theory  to  marine  propellers,  Davidson 
Laboratory  has  sought  to  improve  the  mathematical  model  and  techniques  for 
solving  the  surface  integral  equation,  and  evaluating  the  propeller  loading 
distributions  and  resulting  hydrodynamic  forces  and  moments  (see  for 
example  Refs.  1-5).  A  new  method^  has  been  developed  for  the  inversion 
of  the  downwash  integral  equation  which  employs  the  "generalized  lift- 
operator"  technique.  This  new  approach  eliminates  many  of  the  numerical 
difficulties  arising  In  the  commonly  used  "mode-collocation"  method. 

Certain  mathematical  simplifications  have  been  introduced  in  these 
studies,  in  addition  to  those  imposed  by  the  linearized  version  of  the 

problem.  The  major  simplification  Is  the  approximation  of  the  hrlicoida) 

2 

propeller  wake  by  a  "staircase"  function.  In  Reference  I,  which  presents 
the  exact  treatment  of  the  heiicoidal  wake  for  a  simple  sector-form  blade 
with  flat-plate  chordwlse  loading  distribution,  the  results  are  found  to 
agree  well  with  those  for  the  same  propeller  computed  with  the  staircase 
approximation  of  the  wake.  Admittedly,  only  one  set  of  calculations  was 
performed,  under  several  conditions;  the  complications  in  the  exact  wake 
treatment  considerably  lengthened  the  computation  time  on  a  high-speed  com¬ 
puter  (IBM7090)  and  appeared  too  formidable  to  allow  utilizing  that  method 
more  fully  at  that  time,  it  was  felt  that  use  of  the  exact  treatment 
should  be  deferred  until  a  time  when  more  accurate  experimental  techniques 
would  be  established  for  measuring  the  hull  wake  upon  which  all  calcula¬ 
tions  depend,  thus  Justifying  greater  accuracy  in  the  theoretical  treatment. 

In  the  meantime,  the  numerical  procedure  In  the  solution  of  the 
integral  equation  was  being  Improved,  a)  by  incraaslng  the  number  of 
chordvlse  modes,  b)  by  introducing  more  appropriate  numerical  schemes  to 
deal  with  the  singular  behavior  of  the  kernel  function,  with  the  problem 
of  truncation  of  infinite  series,  and  with  the  problem  of  "overlapping" 
of  the  wake  generated  by  wide  blades,  and  c)  by  accommodating  arbitrary 
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blade  camber  variation.  In  addition,  the  problem  of  Instability  of  the 
chordwlse  loading  distribution  was  treated  and  stability  achieved.^ 

The  results  of  a  number  of  theoretical  calculations  over  a  wide 

range  of  expanded  area  ratio  have  shown  that  the  trends  of  the  hydrodynamic 

4 

forces  and  moments  follow  the  experimental  trends.  The  magnitude  of 
steady-state  thrust  and  torque  have  been  shown  to  agree  within  5-8$  and 
the  magnitudes  of  unsteady  thrust  and  torque  within  5-15$  with  those  of 
NSkDC  experiments  of  a  series  of  3-blade  Propellers  operating  In  an 
accurately  measured  screen-generated  wake. 

During  the  period  of  these  Investigations  also,  the  computer  program 
has  been  converted  for  the  CDC  6600  high-speed  digital  computer  so  that 
the  running  time  has  been  reduced  considerably.  Now  information  about 
steady-state  and  time-dependent  propeller  loading  distributions  and  corres¬ 
ponding  hydrodynamic  forces  and  moments  about  blade  bending  moments  and 
pressure  signatures  at  various  field  points  can  be  obtained  in  a  compara¬ 
tively  short  time. 

Because  of  the  availability  of  high  speed  digital  computers  the  task 
previously  considered  formidable  has  been  undertaken,  viz.,  to  Improve  the 
mathematical  model  by  removing  the  so-called  "staircase"  approximation. 

A  new  kernel  has  been  developed  which  takes  into  account  the  exact  he II— 
coidal  wake  and  the  exact  blade  geometry  through  the  entire  range  of 
expanded  area  ratio.  Thus  an  exact  linearized  version  of  the  unsteady 
lifting  surface  theory  as  adapted  to  the  propeller  case  is  established. 

A  numerical  procedure  and  corresponding  computer  program  is  developed 
here  which  will  furnish  the  blade  loading  distribution  and  the  resulting 
hydrodynamic  forces  and  moments  at  any  desired  frequency.  In  addition, 
the  program  has  the  capability  of  evaluating  blade  bending  moments  about 
the  pitch  line  at  any  radius  and  the  pressure  field  around  the  operating 
propeller. 

This  study  was  sponsored  by  tho  Naval  Ship  Research  and  Development 
Center,  Hydromechanics  Research  Program,  Contract  NOOOl4-67-A“"Q202-OOl8. 
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LINEARIZED  UNSTEADY  LIFTING-SURFACE  THEORY 

The  linearized  unsteady  lifting-surface  theory  for  a  marine  pro¬ 
peller,  with  Its  blades  lying  on  a  hellcoldal  surface  and  operating  In 
nonuniform  flow  of  an  Incompressible,  Ideal  fluid,  was  formulated  by 
means  of  the  acceleration-potential  method  In  previous  Davidson  Laboratory 
papers.  The  development  presented  here  Is  an  adaptation  of  the  previous 
developments,  treating  the  hellcoldal  wake  exactly  as  In  Ref.  1,  and 
incorporating  all  the  Improvements  of  the  later  papers,  the  exact  blade 
geometry,  a  large  number  of  chordwlse  modes,  the  "generalized  lift-operator" 
technique. 

It  Is  known  that  evaluation  of  the  blade  loading  distribution  requires 
the  solution  of  the  following  surface  integral  equation  which  relates  the 
known  velocity  distribution  to  the  unknown  blade  loading, 

W(x,r,cpjt)  JJ  L*(?,p,0;t)  K(x,r,cp;5,p,ejt)dS  (0 

S 

where  x,  r,  9  and  ?,  p,  9:  cylindrical  coordinates  of  control 

and  loading  points,  respectively 

t:  time,  sec. 

2 

S:  propeller  surface,  ft 

W:  known  velocity  distribution  normal 
to  propeller  blades,  ft/sec 

2 

L* :  unknown  loading  distribution,  lb/ft 

K:  kernel  function  representing  the  self- 
induced  velocity  at  a  point  on  the  blade 
due  to  unit  amplitude  lead  at  each  and 
every  point  on  the  blade,  ft/lb-sec 

The  oscillatory  loading  and  velocity  distributions  can  be  written 
as: 
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L'(S>p,ejt)  *  2  L«(x>($,p,e)e,xnt 
\c0 

(2) 

W(x,r,q>jt)  =  £  v/q)(x,r,<p)elqnt 

q=0 

(3) 

where  X  and  q  are  positive  integers  and  designate  order  of  harmonic  of 
the  loading  distribution  and  of  the  inflow  field,  respectively,  the  former 
to  be  determined  by  the  analysis.  The  propeller  rotates  with  angular 

*d» 

velocity  -ft.  The  velocity  W  Is  caused  by  flow  disturbances  such  as  those 
due  to  hull  wake,  Incidence  angle  or  camber  of  the  blades.  In  the  linear 
theory  their  various  effects  are  treated  separately  and  simply  added 
together. 

The  kernel  function  K  has  been  developed  in  Ref.  6  without  resorting 
to  any  geometric  approximation.  In  that  reference,  which  Is  concerned 
with  two  interacting  propellers  of  counterrotating  or  tandem  propulsive 
systems,  a  numerical  procedure  was  suggested  but  not  carried  through  to 
a  computer  program.  Furthermore  treatment  of  numerical  difficulties, 
such  as  the  Cauchy-type  singularity  in  the  separable  form  of  the  kernel 
and  a  higher-order  singularity  with  finite  "Hadamard"  contribution  In  the 
spanwlse  integration,  was  merely  outlined  or  referenced.  The  present 
development  follows  that  of  Ref.  6  but  with  greater  attention  to  detail. 


DEVELOPMENT  OF  THE  KERNEL  FUNCTION 


It  is  known  that  the  pressure  field  generated  by  a  lifting  surface 
S  is  given  by  distributed  doublets  with  axis  parallel  to  the  local  normal, 
and  with  strength  equal  to  the  pressure  Jump  across  the  surface  S.  Thus 
the  pressure  P  at  a  point  (x, r,cp)  at  tiiae  t  wil!  be  given  by 


^(x,r,cpjt)  = -i-  JJ  Ap(S,p,9;t)  - - - - - 

^Pf  S  dn  R'(x,r,9;  S,P,8) 


dS 


(*) 


% 


minus  sign  for  right-hand  rotation  (0  has  absolute  value) 


k 
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In  Eq.  (4),  Y  Is  the  acceleration  potential  function  defined 

as  P/pf 

Pj.  fluid  mass  density 

3/3n  normal  derivative  on  the  surface  S  at  the 

loading  point  (|,p,0) 

•4 

n  unit  normal  vector  having  positive  axial 

component 

Ap(5,p,8;t)  pressure  jump  across  the  lifting  surface, 

l.e.,  Ap  =  P  -  P  ,  pressure  difference 

T  “ 

between  positive  and  negative  oriented 
surfaces  (with  respect  to  x-axts) 

R'(x,r,cpj§,P,e)  C(x-?)2  +  r2  +  pa  ~  2rpcos(9-cp)];- 

e  Descartes  distance  between  the  given 
control  point  and  loading  point 

For  doublets  with  pulsating  strength  Ap(§, p,0)e*U)t  at  point  (5,p,8) 
which  rotates  with( angular  velocity  -n>  Eq.  (4)  yields 

=  J--  JI  - ! - «  (5) 

pf  s  *•(*,  r,<pj5,p,eo-nt) 

where  10  »  frequency,  radlans/sec 

9q  s  angular  position  with  respect  to  the  blade  midchord  line 

When  the  relation  between  velocity  potential  $  and  acceleration 
potential  Y  Is  utilized,  viz., 

*(x,r,cpjt)  «  -  -  |  r,cp; t  -  -2jp^)dT«  (6) 

"CO 

where  U  «=  forward  velocity  and  the  lifting  surface  is  identified  as  t'  i 
helicoidal  surface  of  an  N-biaded  propeller,  where  both  control  and  loading 
points  rotate  with  angular  velocity  -Q  ,  then  the  expression  for  the 
velocity  potential  Is  given  by 
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*(x,r,<P0  ;t) 


n  *  jxnt 
H  2  2 — 

n«l  x-0  wpfU 


SS  AP(X)(5,P,0o)  Jw  elXCa(T»-x)-Sn]  (^)dT,dS  (7) 


where  a  *»  fl/U 

X  *=  u>/fl  ,  order  of  harmonic 
0r  «  2tT/N  (n-1),  n  «=  1,2. ...N 

R  «=  i(',,-5)a  +  r2  +  p2  -  2rpcos[6o-{po+9n-a(T»-x)3}^ 

<p  *  angular  position  of  control  point  with  respect  to  blade 
mid-chord  line 

The  self-induced  velocity  at  (x,r,cf>0;t)  on  the  hellcoidal  surface 
(t.e.,  R.H.  side  of  Eq.  (I)),  will  be  given  by 


b 


N  <*>  HOt 
ID  S  e  * 

n=l  teo  WfU 


s 


a)(5,p,e, 


\  d  I*  IX[a(r»-x)-9  3  3  /  /q\ 

»>5n*’  J  •  "  ¥TdT  dS  (8> 


where  Is  the  normal  derivative  on  the  hellcoidal  surface  at  (x,r,tp), 
the  control  point. 

The  directional  derivatives  normal  to  the  hellcoidal  surface,  which 
is  given  by  x  «=  cp^a  or  §  «  0^/a  ,  are: 


3  r  /_  9  1  &  \ 


(9) 
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(!o: 

The  L.H.  of  Eq.  (l),  viz., 

»  .  .  ,  _  »  .  N  |q(Clt-«p  ) 

W(x,r,cp;t)  >  S  W<q>  (x,r,<P  )e,q  *  =  £  V^q\r)e  °  (11] 

°  q=0  °  q=0 

which  Is  the  known  velocity  distribution  obtained  from  measurements,  has 
time  dependence  of  the  form  e*q^  ,  whereis  the  right-hand  side, 
which  Is  equivalent  to  Eq.  (8),  has  time  dependence  of  the  form  e^*. 
Thts  dictates  that 


X.  =  q 


hence 


-  J*  J  L,(q)  (P,0o)e,qClt  K(r,<po;p,eo;q)dS 


where 


L‘(q)(p,9j  s  AP(X)  (p,9j,  lb/ft2 


-  ■*{?  i,  -'"5"  ^  &  <>■ 

CD  -  © 

8  «  -2-.-SL  -  (x  -  O 


and  by  8  -»  0  Is  meant  x  -*  «pQ/a  and  5  "*  8  /*•  The  limiting  process  Is 
Introduced  to  avoid  the  difficulty  arisinr  In  the  mathematical  manipulation 
due  to  the  hlgh-order  singularity. 

The  surface  integral  given  by  Eq.  (12)  Is  equivalent  to 


7 


R-1509 


where  0^  Is  projected  semichord  length  of  the  propeller  blade  at  the  load¬ 
ing  point,  In  radians,  and  the  factor  (Vi  +  a2pa/ap)  Is  the  result  of 
changing  the  Integration  over  the  actual  propeller  blade  to  Integration 
over  Its  projection  In  the  propeller  plane.  It  should  be  noted  that  this 
does  not  correspond  to  placing  the  dipoles  on  the  projection  of  the  blade 
in  the  x  a  0  plane. 

With  the  transformation 

eo  =  aP  -  eg  cos0a  (15) 

whwere  0^  is  angular  chordwise  location  of  the  loading  point,  and  Is 
angular  position  of  the  midchord  line  of  the  projected  blade  from  the 
generator  line  through  the  hub  (<7P  is  a  measure  of  skewness  and  Is  positive 
towards  the  trailing  edge)  and  letting 

l(,)  (p,e0)  *  L'(q)(p,e0)  peg  (it/ft)  (u>) 

the  Integral  becomes 

1  ■  L  J‘{'-(q)( p.9a)*,<,ntK}  slnB^  dp  ( 17) 

Op  p 

Use  will  be  made  of  the  following  expansion  scheme  for  Hie  inverse 
Descartes  distance 

«=  l/{(T,-$)2  +  r2  +  pz  -  2rp  cose}5 

S  J*  lm(lklp)  Km(lklr)e,(T,"^k  dk  (18) 

m=  - 00  -00 
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for  p  <  r  ,  where  |  (  )  and  K  (  )  are  modified  Bessel  functions  of  the 

ni  m 

first  and  second  kinds.  For  p  >  r,  p  and  r  are  Interchanged  In  the 

product  of  the  modified  Bessel  functions.  Here  3=6  -  cp  +  5  -  a(T*-x). 

o  o  n 

After  taking  the  derivatives  as  indicated  In  Eq.  (13)  and  Integrating 
over  Ti  (see  Ref.  I),  noting  that 

£  rN  for  (n>-q)  -  M  Uo,  ±1,  «... 

"  6  n  =  ip  for  all  other  m-q 

n=1 


the  kernel  function  becomes  for  p  <  r 


K(r#9o,P,0ojq)  »  -  tjttTTJ 


pr 


Hm  £  tm(9o-?o) 

Pfu  [(l+aap2)(l+a2r2)]^  ^  m-  -•  # 

m=q+  JtN 


.  {.,a(”-<l)(K-5>C.2(ni-q)  +  4  ]  Ca2(m-q)  +  -—]  lm(alm-qlp)  Km(alm»qlr) 


“  *  J  <ak  +  7^ak  +  *3)® 

-09  r  P 


r2'  r 

lk(x-5)  "JIMP 

k  -  a(m-q) 


<fk  }  ( 19) 


It  Is  obvious  from  Eq.  (13)  that  In  the  limit  as  6  -*  0,  I.e.,  as 

a  -*  ax  and  0  -*  a5,  R  will  go  to  zero  when  n  B  r,  3  =0  and  T*  =  0. 
tq  o  r  n 

Appendix  A  Investigates  this  singular  behavior  of  the  kernel  before  and 
after  the  l/R  expansion  of  Eq.  (18)  and  shows  that  the  singularity  Is  of 
high  order  with  "Hadamard"  finite  contribution. 

In  addition,  a  Cauchy-type  singularity  is  introduced  in  the  k-integra- 
tlon  of  Eq.  ( 19)  through  the  series  expansion  of  l/R,  which  nevertheless 
is  essential . for  the  reduction  of  the  surface  integral  equation  (IT)  to 
a  line  integral  equation.  The  kernel  function  is  now  in  separable  form 
which  facilities  the  chordwise  Integration  and  permits  use  of  the 
"generalized"  lift  operator  technique  for  Inverting  the  Integral  equation. 
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REDUCTION  TO  A  ONE-DIMENSIONAL  INTEGRAL  EQUATION 

The  unknown  loading  L^(p,8a)  of  Eq.  (IT)  will  be  approximated  In 
the  chordwlse  direction  by  the  BIrnbaum  mode  shapes 

L(q)<p,6j  =  i  {l(q'°(p)  8(1)  +  S  L(q*S)  (p)  8(5)}  (20) 

n»2 

where  L^#"\p)  are  the  spanwlse  loading  components  to  be  determined  by 
the  solution  of  the  integral  equation  and 

e 

0(  1 )  e  COt  -jp 

@{n)  *  sin  (n-i)8a  ,  n  >  i  (2i) 

Now  after  the  limit  is  taken  in  Eq.  ( 19)  and  with  the  trigonometric 
transformation  of  (15)>  the  9^- integration  can  be  performed  analytical 1y. 

Both  sides  of  the  downwash  integral  equation  are  then  operated  on 
by  the  "generalized”  lift  operator,  the  structure  of  which  Is  dictated  by 
the  separable  form  of  the  kernel  (see  Refs.  2  and  5).  This  method  elimin¬ 
ates  the  cpQ-dependence  of  the  integral  equation.  Thus  a  set  of  line  integral 
equations  is  obtained  with  maximum  order  of  lift  operator  m  equal  to  maxi¬ 
mum  number  of  unknown  chordwise  modes  n.  The  solution  of  the  m  =  n  Integral 
equations  is  obtained  by  the  collocation  method. 

With  the  transformation 

*0  ■  qr  -  ®b  ""fa 

where  tp  is  angular  chordwlse  location  of  the  control  point  and  the  super¬ 
script  r  refers  to  values  at  the  radial  location  of  the  control  point, 
the  non-dimensional  I zed  downwash  integral  equation  is  reduced  to  the 

following  form  for  each  q  ^  0  and  each  m  =  1,  ...  n  : 

3  n  ’  max 
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f  -  N  \  _r  yv  -ImAo 

*■  4ttp. r  avi+a^  m=  -  « 
f  i^q+JtN 


.  {•»(««l)to[a8(»-q)  +  -|3  [a2(m-q)  + -%]  l(5)(qeJ)A(5)(qeJ) 

r  P 

.  I  (alm-qlp)  K  (alm-ql r) 
m  m 


-  ^  J  (ak+  4y)(ak+  -lyje1  a^ 


l(S)((n> 


,  X>A 


(n 


~)e£) 


k-a(m-q) 


lm(lklp)Km(lklr)dk} 


(27) 


/-\  .  TT  -IxCOSB 

where  A^n'(x)  “  ^  J  0(n)e  slnB^dB^  which  is  also  evaluated  in 

o 

Appendix  B.  In  Eq.  (27)  all  terms  outside  the  first  braces  are  non- 
dimensionalized  with  respect  to  rQ,  as  is  also  p  in  Eq.  (22),  and 
Ao  *=  <jP  - 

Let  X  »  k  -  aj&N,  then  the  k-lntegral  of  Eq.  (27)  Is 


jjIj&nao  p®  g(x)dx  m  e  unao  p®  aMraBl  d\ 

*  .1  x  x 


(28) 


where 


1  4 


Act 


•<*>  "  <IK>m  Bm,n<X>e 

|>m(IX+a£Nlp)  Km(IX+ajeNlr)  for  p  <  r 


OK) 


m 


^m(IX+aXNIr)  Km( I X+a^N I p)  for  p  >  r 


Bm,n (x)*  (aX+a2jeiH- 


2,)  (ax+.’iis  +  4)i(S)((q-  i)eJ)A(5)(h-  i)eg) 

«*  P 
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Equation  (28)  has  an  integrable  singularity  at  \  ■  0.  The  value  of  the 
integrand  at  this  singularity  is  obtained  by  (.'Hospital's  rule,  as  shown 
in  Appendix  C: 

a)  when  \  ■  0  and  i  ¥  0 

IS  * 2  (<«>.,  J*  v«»  *  ^  I J 


a  Ok). 

i  B-  -(0)  -—2 

m,n'  o\ 


V 


(29a) 


where  the  upper  sign  is  taken  when  &  >  0  and  the  iowe.  ^hen  l  <  0,  and 
the  derivatives  are  given  in  Appendix  C. 


b)  when  \  *  0,  l  *  0  and  m  =  q  *  0 


I  im  * 
X-0  J 


q  -  o 


(29b) 


c)  when  \  »  0,  1*0  and  m  «  q  +  0 


L]  -  *  (IS  (.«.!  J 


where 


r  o  *  P 


im4 

’  177 


[sji 1 <i5) (<je£)A<?) (qej)-*jl (S) (qerb)A{f') (qej)]}  (29c) 


ltai(IK)  | 

tr*0  mlx-0 


U^'"' 

TOT  ©'■'  *«■  »  >  ' 


and  I  (x)  and  Aj^(x)  are  defined  in  Appendix  B. 
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Equation  (27)  can  now  be  written  as 
.(m»n)„  / . -_N 

rr  .  J 

m-q+jtN 

where 


{_^_}  L±^  £  {9(0,  -i  f  (30 

l4tTpfU2r  J  a/l+a2r2  "i— ® 

»  0  MMiix  «u 


g(0)  ■  (IK)  I  { (a2£Ht  ^)  (a2AN+ *^)  I  ^  (qO^M^  (qB^)} 

X-0  u  r  p 

For  large  I x I  *  Up,  the  integrand  of  Eq.  (30)  tends  to  zero. 

For  large  argument  l±X+ajjNl(p  or  r)  of  the  modified  Bessel  functions, 

-liA+aiNI  lr-pl 

(IK)  of  (28/  *  £_ - 

2  ,/r£  l±X+aANi 


so  tha* 

g(^)  e 


l±X+a£NI  lr-pl  2  .  m  _ _  m  ^±i  ho 


a4(±X+ajtN  +  -^r)  (±X+aAN+ ~2*)e  a 
*•  2/rp  xiiA+aiNI  ar  ap 

.  l<l5)((qii)9D  A(f!)((q?M)  <3') 


If  x  and  a£N  are  opposite  in  sign  and  iXHIa^NI,  i.e.,  the  arguments  are 
sma 1 1 ,  then 

1  Imt 

^m5*  21ml  1 

where  z  *  p/r  for  p  <  r  and  z  «=  r/p  for  p>r.  In  this  case 


tfcaU-L. 

X  21ml 


zlml 


(i)((q+"  i)0^)A(")((<f3  ileCi  (32) 


“TT 

xr  P 


When  x  is  large,  the  product  of  the  l^and  A ^ functions  varies  as  1/X. 
Equation  (31)  is  seen  to  be  negligibly  small  for  large  \  provided  that 
Iml  -  Iq+j&NI  is  not  too  much  larger  than  for  example  if  axN  is  smaller 
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2 

than  x  •  Equation  (32)  tends  to  zero  for  both  large  \  and  large 
iml  £  Mp  since  I \l  -*  lajtNI. 

The  Integral  equation  given  by  Eq.  (22)  with  the  kernel  given  by 
Eq.  (32)  can  be  solved  numerically  by  the  usual  collocation  method.  In 
this  method  the  loading  L^,n^(p)  Is  assumed  to  be  constant  over  each 
small  radial  strip,  23  In  length.  Then  only  the  kernel  needs  to  be 
integrated  over  the  radial  strip.  When  p  /  r  the  integration  of  the  kernel 
In  the  range  pj  -j3<p<Pj  +  pls  performed  by  the  tangential  method. 

In  the  radial  strip  r-p<p<r+p  which  Includes  p  **  r,  a 

high-order  singularity  occurs,  as  noted  in  Appendix  A.  Its  "Hadamard" 

finite  contribution  Is  evaluated  by  the  Lagrange  Interpolation  method 

c 

described  in  Appendix  D,  as  suggested  by  Mangier. 

PROPELLER  LOADING  AND  RESULTING  FORCES  AND  MOMENTS 


Loading 

_  The  downwash  integral  equation  (22)  is  solved  for  propeller  loadings 
L *n)  (p)  due  to  the  known  velocity  distribution  W^>m^(r)  normal  to  the 
blades.  The  velocity  distribution  of  Eq.  (26)  is  the  hull-induced  dis¬ 
tribution  at  shaft  frequencies  q  obtained  from  wake  surveys  in  the 
propeller  plane.  To  the  loadings  due  to  thesn  wake  velocitlc.  are  added 
loadings  arising  from  various  disturbances  Imposed  on  this  flow. 

The  propeller  blades  do  not  coincide  with  the  assumed  helicoidal 
surface  of  pitch  i/a  but  are  located  on  a  nearby  surface  of  pitch  P(r) 
and  furthermore  they  are  cambered.  (As  in  previous  Davidson  Laboratory 
papers,  thickness  is  neglected.)  The  velocities  induced  by  the  flow- 
incidence  angle  and  camber  effects  are  independent  of  time  because  the 
blades  are  considered  rigid,  so  that  only  the  steady-state  loading  (q  “  0) 
will  be  affected.  For  the  flow  angle  effect  the  left-hand  side  of  the 
integral  equation  is  replaced  by 

(r)  - - 

- -  -  J\  +  a2r2(tan“l  -  tan"1  ~r)  for  m  »  1,2 

=  0  for  m  >  2 
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(since  this  effect  depends  only  on  radial  position  and  not  on  cpa)  where 

the  incident  flow  angle  is  defined  by  the  difference  between  the  actual 

geometric  pitch  angle  9  (r)  *  tan"' (P/2ttO  and  the  hydrodynamic  pitch 
-  J  P 

angle  tan  (1/ar)  of  the  assumed  surface. 

For  the  camber  effect  the  left-hand  side  of  the  Integral  equation 


c 


JhiiL  r"  $  (m)  d£Xusl 

2rrb(r)  o 


/l*a2r2  d% 

isur  j.  *(b)  se.st«; 


where  b(r)  =  semichord,  ft. 

e  slope  of  camberline  f(r,s)  given  at  discrete  points 
°  measured  from  face  pitch  line 

s  ■  (l-co3cpa)/2,  chordwise  location  nondimensionalized 

on  basis  of  2b (r) 

The  evaluation  of  integral  (34)  Is  given  in  Reference  4  for  arbitrary 
blade  camber  shape. 

Once  the  values  of  L^'n^  due  to  the  various  disturbances  are 
determined  the  spanwise  loading  distribution  follows  from  Eq.  (20): 

l.(q)(r)  -  L^(r,ejsineadett 

"nl  U*q,1V)(l+cose  )  +  £  L(q,n)(r)sin(n-l)o  slnelde 

O  a  pm2  &  Cr  Ot 

-  U(q>l) (r)  +  i  L(q’2)(r)  (35) 

Propeller-Generated  Forces  and  Moments 

With  the  principal  components  of  the  propeller-induced  forces  and 
moments  defined  as  follows 


Forces :  F. 


thrust  (x-direction) 


F  and  Fz  «  horizontal  and  vertical  components, 
Y  respectively,  of  the  bearing  forces 
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Moments:  *  torque  about  the  x-axis 

Q  and  Q  *  bending  moments  about  the  y-  and  z-axis, 
y  z  respectively 

the  total  forces  at  frequency  /N(jfc*0, 1,2,...)  acting  on  an  N-bladed  pro¬ 
peller  will  be  given  by  (cf.  Ref.  3) 

Fx  *  Re  |NrQ  e,,iNnt  (r)cosep(r)drj 

Fy  •  Re  [L(",-,){r)+L<^,)(r)}ln0p(r)dr} 

Fj  .  Re  {£°  .'*»  jl  [L("!-')(r).t<^,>(r)]sln9I>(r)dr} 

(with  the  sign  convention  adopted  In  the  present  investigation). 

The  moments  are  determined  by  * 

Qx  -  -  Re  e,ANnt  Jl  (r)sin6p(r)rdr \ 

Nr2  ° 

Qy  -  Re  fl°-  e1*™  (r)+L<^> (r)]cosep(r)rdr} 

2  0 
Nr 

d2  -  Re  {-^e,1Wit  J’,[L(/N',)(r)-L(j!N+,)(r)]cos0p(r)rdr} 

5 

It  may  be  obser  fed  from  the  foregoing  that  the  propeller-generated 
transverse  forces  and  bending  moments  are  evaluated  from  propeller  load-  j 

Ings  associated  with  wake  harmonics  ar.  frequencies  adjacent  to  blade  3 

frequency,  i.e.,  at  q  *=  ju'il,  whereas  the  thrust  and  torque  are  deter- 

% 

mined  by  the  loading  at  blade  frequency.  The  steady-state  thrust  and  | 

torque  are  determined  at  zero  frequency;  the  corresponding  mean  transverse  5 

forces  and  bending  moments  are  determined  at  shaft  frequency.  | 

Blade  Bending  Moments 

The  blade  bending  moment  about  the  face  pitch  line  at  any  radius 
rj  of  a  blade  is  calculated  from  L^(r)  at  any  shaft  frequency  q  as  * 

H(<!)  .  r2  elqn«  j1  L(q)(r)coj[»p(r)-ep(r.)]  (r-r.)dr  (38)  ;  j 

r  •  I 

I  | 

j  I 

i  5 

*  r 

Z  i 
»  © 

%  | 
i  ? 

M  M 
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The  total  vibratory  blade  bending  moment  distribution  Is  then  obtained  from 

Mb  =  Re  2  Mb(q)  elqnt  (39) 

q 

The  computer  program,  which  implements  the  numerical  procedure 
developed  here,  furnishes  the  blade  loading  distributions  (Eq.  35),  the 
hydrodynamic  forces  and  moments  (Eqs.  36-37),  and  the  blade  bending  moment 
distribution  (Eqs.  38  and  39). 

The  program  also  Implements  the  procedure  developed  in  Ref.  10  for 
the  evaluation  of  the  p  opel ler-induced  vibratory  pressure  field  from  the 
blade  loadings  L^,n^(p).  The  blade-frequency  pressures  due  to  loading, 
are  shown  to  be  of  the  form 


(£N) 


2 


2 

(J,*'0 


2 

p 


{!>•"■)  (p)k(^,!N'")+  conj.  bu(M..n)  (p)kfg.+XN,n)-j j.  {40) 


where  the  k  functions  are  the  so-called  Influence  functions  connecting  the 
loading  points  on  the  blade  with  the  space  points. 


Also  available  from  Ref.  10  Is  a  computational  procedure  to  deter¬ 
mine  the  blade  thickness  effects  on  the  pressure  field.  Although  these 
effects  are  negligibly  small  in  the  case  of  bearing  forces,  they  are 
demonstrably  large  in  the  pressure  field  in  the  neighborhood  of  the 
operating  propel ler. 


NUMERICAL  RESULTS 

Using  the  present  development  with  exact  treatment  of  the  he  II  col  dal 

wake  from  the  blades,  systematic  calculations  of  steady  and  unsteady  blade 

loading  and  hydrodynamic  forces  and  moments  were  made  for  a  series  of 

3-bladed  propellers  of  expanded  area  ratio  EAR  between  0.3  and  1.2.  These 

3  4 

propellers  were  chosen  because  they  were  subjects  of  earlier  papers  ’ 
which  had  used  the  "staircase"  approximation  of  the  wake,  and  because 
experimental  data  for  the  three  propellers  of  EAR,  0.3,  0.6  and  1.2,  were 
available  from  controlled  NSROC  tests  In  nonuniform  flow  with  wake  structure 
rigidly  specified  and  vibratory  forces  made  large  to  minimize  experimental 
error. 
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The  tests  had  been  conducted  In  the  NSRDC  24-Inch  water  tunnel  using 
the  closed-jet  test  section  with  a  screen  to  produce  the  wake,  a  3-cycle 
screen  for  the  force  along  and  moment  about  the  longitudinal  axis  and  a 
4-cycle  screen  for  the  forces  along  and  moments  about  the  transverse  and 
vertical  axes.  The  tests  were  run  at  a  constant  speed  of  15  rps  with 
free  stream  velocity  close  to  12.5  ft/sec.  The  design  mean  thrust  coeffi¬ 
cient  K^.  was  0.150  (practically  open-water  value  In  the  3-cycle  screen 
wake) . 

The  propellers,  with  destroyer-type  blade  outlines  and  modified 
NACA-66  section  with  a  *  0.8  mean  line,  were  one  foot  In  diameter  with 
hub  diameter  0.2  ft  and  pitch  ratio  1.08  at  .7  radius.  For  the  propellers 
at  other  expanded  area  ratios  the  pitch  variation  along  the  span,  RPM  and 
flow  velocity  normal  to  the  blade  were  considered  to  be  the  same  and  equal 
to  the  corresponding  values  for  the  propeller  of  EAR  “  0.6  (taken  as  the 
norm).  The  camber! ine  variations  were  determined  by  linear  Interpolation 
of  the  known  ratios  of  maximum  camber  to  chord  m  /c  shown  in  the  table 
below  for  the  propellers  with  cAR  •  0.3,  0.6  and  1.2. 

TABLE  1 

Ratio  of  Maximum  Camber  to  Chord,  mx/c 


Radius 

EAR  «  0.3 

EAR  »  0,6 

EAR  »  1,2 

.25 

.0370 

.0228 

.0195 

.35 

.0388 

.0231 

.0202 

.45 

.0372 

,0224 

.0196 

.55 

.0340 

.0212 

,0185 

.65 

.0312 

.0203 

.0177 

.75 

,0290 

.0198 

.0170 

.85 

.0270 

.0189 

.0160 

.95 

.0247 

.0174 

.0147 

The  results  of  the  calculations  are  shown  In  Figures  1-4  and 
Tables  2  and  3.  Figure  I  compares  the  mean  thrust  coefficients  evaluated 
by  the  present  method  with  the  earlier  computations  using  the  "staircase" 
approximation  of  the  blade  wake,  and  with  the  experimental  value  of  0,15. 

As  can  be  seen,  with  the  exception  of  the  propeller  of  0,3  EAR,  the  present 
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calculations  are  closer  to  the  measurements.  Repeating  the  calculations 
with  finer  mesh  and  higher  limits  for  the  integration  has  not  changed  the 
result.  The  approximately  13%  increase  in  mean  thrust  at  EAR  *  0.3  re¬ 
flects  the  disproportionately  large  camber  employed  with  the  small  blade 
area  ratio  (see  Table  1). 

It  is  to  be  remembered  that  the  blade  thickness  effect  on  the  loading 
has  been  neglected  here  as  in  the  previous  studies.  This  effect,  which 
Is  expected  to  decrease  the  calculated  mean  thrust,  is  probably  insigni¬ 
ficant  in  most  cases,  but  may  be  large  enough  in  the  case  of  the  unusually 
thick  propeller  of  0.3  EAR  to  bring  the  value  closer  to  experiment.  This 
is  the  subject  of  a  current  investigation  at  the  Davidson  Laboratory. 

The  steady-state  spanwise  loading  distributions  are  presented  in 
Figure  2  in  terms  of  EAR.  At  all  EAR  the  loading  is  shown  to  increase 
with  radius  from  0.25  to  0.85  and  then  decrease  at  0.95.  The  trends 
versus  EAR  at  every  radius  with  the  exception  of  0.95  aro  consistent. 

The  rise  in  all  curves  at  EAR  ■  0.3  is  reflected  in  the  higher  mean 
thrust  coefficient  for  that  propeller. 

Figure  3  compares  the  blade-frequency  vibratory  thrust  and  torque 
amplitudes  computed  by  the  present  method  with  the  earlier  calculations 
reported  in  Reference  4  which  used  the  "staircase"  approximation  of  the 
blade  wake,  and  with  the  experimental  values.  The  present  theoretical 
results  are  quite  close  to  the  test  data.  The  proportionately  greater 
improvement  in  the  torque  amplitudes  over  the  earlier  calculations  is 
especially  gratifying  since  it  Indicates  that  the  radial  load  distributions 
are  more  accurately  determined. 

Figure  4  presents  the  real  and  imaginary  parts  of  these  redial  load¬ 
ing  dlstribut Ions  versus  EAR  at  first  blade  frequency  q  *  N  «  3.  The 
distributions  show  smooth  trend*  with  EAR,  and  a  definite  pattern  with 
increasing  radius  from  0.25  to  0.85.  At  0.95  radius  the  trend  is  completely 
different.  The  curves  of  loading  es  a  function  of  EAR  are  similar  to 
those  shown  In  Ref.  4,  Figure  10,  but  are  lower  in  magnitude  and  without 
the  inflections  at  high  EAR  of  the  earlier  calculations. 
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A  comparison  of  the  blade-frequency  amplitudes  of  nondimensionai 
transverse  and  vertical  force  coefficients  and  coefficients  of  moments 
about  the  transverse  and  vertical  axes  from  NSRDC  measurements  and  the 
theory  is  given  in  Table  2.  The  coefficients  are  for  the  3-biaded  pro¬ 
pellers  tested  in  a  4-cycle  screen-generated  wake  at  the  design  advance 
ratio  J.  (The  resolution  of  forces  and  moments  Is  shown  in  Fig.  5.) 

TABLE  2 

Experimental  and  Theoretical  Transverse  and  Vertical  Force  and  Moment 
Coefficients  for  a  Set  of  3-Bladed  Propellers 


EAR 

-  0.3 

EAR 

■  0.6 

EAR 

-  1.2 

2 

Amp.  x  10 

Exp. 

Theory 

Exp. 

Theory 

Exp. 

Theory 

kf 

y 

1.08 

0.96 

1.34 

0.98 

0.77 

0.40 

# 

*** 

\ 

1.04 

0.95 

1.25 

1  ;06 

0.57 

0.69 

A# 

s 

0.92 

0.87 

1.10 

0.93 

0.77 

0.40 

\ 

0.89 

0.83 

1.06 

0.93 

0.70 

0.53 

A  comparison  of  the  corresponding  phase  angles  is 

given  In 

Table  3 

which  also  presents 

the  phase  angles  for  blade-frequency  thrust  and  torque 

TABLE  3 

cp »  Phase  Angle  in  Degrees 

■ 

EAR 

-  0.3 

EAR 

■  0.6 

EAR 

-  1.2 

Vibratory  Forces 

Exp. 

Theory 

Exp. 

Theory 

Exp. 

Theory 

A* 

Fx 

207 

222 

238 

257 

323 

331 

M 

<*x 

205 

42* 

234 

77* 

310 

151* 

A# 

Fy 

351 

153* 

32 

204* 

164 

332* 

356 

148* 

31 

191* 

100 

288* 

W 

Fx 

255 

255 

293 

308 

77 

67 

A# 

262 

247 

294 

292 

2 

34 

Approximately  180  deg,  out-of-phase  with  experiment. 
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Knowing  the  magnitude  H  and  phase  angle  9  of  a  particular  force  or 
moment  its  cosine  and  sine  components  can  be  determined  by  means  of  the 
expression 

M  cos  (qfit  -  cp) 

This  form  is  aiso  useful  in  establishing  the  lead  or  tag  of  the  quantity 
with  respect  to  the  sinusoidal  input  wake  at  a  specified  radial  position 
or  the  location  of  a  maximum  relative  to  the  rotating  blade. 

For  all  the  transverse  and  vertical  force  and  moment  coefficients 
the  trends  with  EAR  of  measured  and  theoretical  values  are  similar.  How¬ 
ever,  the  correlation  between  experimentally  and  theoretically  derived 
coefficients  is  not  as  good  as  that  for  the  vibratory  thrust  and  torque. 

Examples  of  propeller  blade  chordwise  loading  distributions  are 
shown  in  Figures  6  and  7  for  the  two  extremes  of  expanded  area  ratio, 

0.3  and  t.2.  These  figures  compare  the  distributions  obtained  by  the 
normal  sum  of  the  Blrnbaum  chordwise  modes  (up  to  finite  n  *  M)  viz. 

S-  -  L(,)(r)  cot  V+  S  L(S)(r)  sin(n-l)0  (41) 

n  n-2  01 

with  the  Cesaro  sum 

a.  -  L(,)(r)  cotf-*P  (2^)  L(5)(r)  sln(5-l>0  (42) 

n  1  n-2  M  a 

The  Cesaro  summability  is  a  proper  procedure  for  obtaining  the  limit 
of  slowly  convergent  or  even  divergent  series.  Use  of  this  method  In 
Ref.  4,  where  calculations  of  loading  were  based  on  the  "staircase"  approxi¬ 
mation  of  the  blade  wake,  served  to  stabilize  the  chordwise  distribution 
which  was  particularly  erratic  In  the  neighborhood  of  the  leading  edge 
for  large  EAR  (>1). 

The  Cesaro  and  normal  sums  in  Figs.  6  and  7  are  both  smooth  and  close 
In  value.  This  is  an  indication  that  removal  of  the  "staircase"  approxi¬ 
mation  with  the  present  exact  treatment  of  the  blade  wake  results  In  a 
more  rapidly  converging  series,  and  consequently  In  a  smoother  chordwise 
distribution* 
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CONCLUSION 

The  theory  which  was  developed  in  Ref.  2  in  adapting  the  unsteady  j 

lifting  surface  theory  to  the  marine  propeller  case  has  been  revised  by  j 

removing  the  so-called  "staircase"  approximation  describing  the  blade  wake.  j 

The  present  "exact"  mathematical  model  takes  cognizance  of  the  exact  he  1 i  — 
coidal  geometry  of  the  propeller  blade  and  its  wake,  as  well  as  the  I 

presence  of  interacting  blades  and  the  fact  that  the  propeller  operates 
in  nonuniform  flow  (hull  wake).  A  new  numerical  procedure  has  been  devised  ; 

and  programmed  for  the  CDC  6600  high-speed  digital  computer  primarily  for 
the  calculation  of  the  blade  loading  distributions  at  various  shaft  fre¬ 
quencies.  The  resulting  hydrodynamic  forces  and  moments  are  evaluated 
from  the  loading  distributions:  thrust  and  torque,  steady  and  time- 
dependent,  from  the  loading  at  zero  frequency  and  multiples  of  blade 
frequency,  the  other  bearing  forces  and  moments,  transverse  and  vertical,  * 

from  the  loading  at  first  shaft  frequency  for  the  steady-state  case  and 
from  the  loadings  at  one  above  and  one  below  blade  frequencies  for  the 

i 

time-dependent  cases.  In  addition,  the  program  has  the  capability  of 
calculating  the  vibratory  blade  bending  moment  about  the  face  pitch  line 
at  any  radius  and  the  pressure  field  due  to  the  loading  of  the  propeller 
when  It  operates  In  uniform  or  nonuniform  inflow  fields. 

The  results  of  this  new  approach  with  the  exact  helicoidal  blade 
wake  have  not  shown  any  dramatic  changes  from  the  results  of  previous 
calculations.  However,  it  can  be  stated  that  an  overall  improvement  has 
been  witnessed  In  the  correlation  of  forces  and  moments  with  experimental 
measurements,  which  indicates  that  the  radial  loading  distributions  are 
more  accurately  determined.  The  detailed  loading  distributions,  both 
radial  and  chordwise,  are  much  smoother  than  obtained  in  the  previous 
calculations  by  means  of  the  "staircase"  model. 

Other  unreported  calculations  for  propellers  In  the  wake  of  a  ship, 
rather  than  in  the  strong  blade-frequency  wake  generated  by  a  screen, 
have  shown  little  difference  in  the  results  by  both  methods  except  in 
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the  steady-sUte  case  where  the  new  method  yields  results  closer  to  the 
experimental  values. 

Unfortunately  the  machine  running  time  has  Increased  considerably, 
calculations  taking  about  4  to  5  times  longer.  It  is  felt,  at  this  stage, 
that  numerical  simplifications  and  other  approximations  should  not  be 
attempted  before  extensive  use  of  this  program  can  be  made  for  a  better 
understanding  of  the  importance  of  the  various  contributions  to  the 
solution. 

It  is  suggested  for  the  present  that  if  detailed  and  accurate 
information  is  required,  the  program  of  the  "exact"  mathematical  mode) 
should  be  utilized,  whereas  when  integrated  effects  like  forces  and 
moments  are  desired  the  "staircase"  model  ts  adequate.  This  statement 
applies  especially  when  the  relative  merits  and  disadvantages  of  propeller 
settings  are  under  consideration. 
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FIG.  2  STEADY-STATE  SPANWISE  LOADINGS  FOR  3-BLADED  PROPELLERS 


FIG.3.  AMPLITUDES  OF  VIBRATORY  THRUST  AND  TORQUE  FOR  3 
BLADED  PROPELLERS 


FIG. 7.  CHORDWISE  DISTRIBUTIONS  AT  0.65  RADIUS  FOR  PROPELLER 
OF  EAR -1.2:  COMPARISON  OF  NORMAL  SUM  AND  CESARO  SUM 
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Appendix  A 


Singularity  of  the  Kernel  Function  at  p  ■  r 


I  Before  the  1/R  Expansion 

If  the  substitution  t  ■  t'  -  §  Is  made  in  Equation  (13)>  the 
kerne)  becomes 


(^PrU)  K 


1  im 
6-«0 


N 

£ 

n-1 


-iq  e. 
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a_ 

an 


x-§ 
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iqa(T-x+§) 


dr 


where  R  “ 
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P2  -  2rp 
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CD  + 
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(t-x+5) ]} 
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(A—  1 ) 


In  the  limit  as  6-4)  ,  R  will  go  to  zero  when  p  -  r  ,  6n  ■  0  ,  and 
t  »  0  .  Therefore  the  singular  behavior  occurs  when  n  ■  1  and  x  -  5  >  0. 
For  a  study  of  this  behavior,  only  that  part  of  the  T-Integratlon  need  be 
considered  that  is  within  a  small  range  around  t  **  0  ,  viz.,  -  y  £  x  £  y 
where  y  Is  a  sufficiently  small  but  positive  fixed  number.  Thus  part 
of  the  kernel,  for  n  *  1  ,  may  be  designated  by  Kj 


-  1  Im 
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J  an' 

-Y 


a_ 

an 


Iqa(T-x+5) 


‘1 


_r 


After  taking  the  derivative*  and  limit':  Kj  becomes 
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+  («2P2  j  0  (agrg  1  0  cos  ar  _  3(a2pg  +  1)  (a2r2  4-  1)  slngaT  d1 
rp(r2  +  r2  +  p2  -  2rp  cos  at  )3/2  (t2  *  r2  +  02  -  J?rp  cos  ar)5/2  j 

(A-3) 

Let  p  ■  r  +  c  where  c  Is, very  small.  Then  the  expressions 
Involved  In  Kj  can  be  expanded  In  a  series  of  »  and  r  .  After  lengthy 
algebraic  manipulation  and  trivial  Integrations,  the  singular  behavior  of 
Kj  Is  described  by  the  following  expression: 

K.  «  -  e~,qa<x“5)  _ — _  fg  A-+  a2r2  +  4a2  1 

1  v  a2r2  \r2(p  -  r)2  r  -A  +  a2r2  p  -  r 

(as  p  -»  r) 

The  first  term  has  a  hlgh-order  singularity  with  "Hadamard"  finite 
contribution  which  Is  evaluated  by  the  method  of  Mangier.^  The  second 
term  Is  the  well-known  "Cauchy"  type. 

1 1  After  the  1/R  Expansion 

For  the  case  of  n  ■  I  ,  x-S  >  0  ,  Eq  (19)  becomes  after  the  limits 
have  been  taken 


*  [(l+a2p2)  (!ra2r2)]^2 


where  (IK)m  »  !„,(  j  k  |  p)  (  (k  |  r) 
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When  |  m  |  ^  |  H  |  large,  the  generalized  mean-rvalue  theorem  can  J>e 
used  for  the  Integral  term  and  this  In  conjunction  with  the  residue  theorem 
leads  to 

d  d 

P  f(k)p(k)dk  «f(A)  f  p(k)d  cSAsd 
c  c  p(k)  *  0 


with  p(k) 


elk(x-5) 
k-a (m-q) 


Ik(x-S) 


*  elk(x-?) 


J„,(k>  ~k-?CT  dk  ~  f(=(m-q))  J„  '  k-aft-q)  dk 


(A-6) 

for  (x-?)  >  0 


Thus  for  large  |mj  the  Integral  term  Is  approximately  equal  to 
the  closed  term  and  with  |  m  j  ^  )  M  J  »q  Kj  reduces  to 


K  _ _ ££ _  Z  -laq(x-S) 

1  [(l«V)(lMar2)3,/2 
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n 

Nicholson'  showed  that  the  product  lm(amp)Km(amr)  can  be  approxi¬ 
mated  by 
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for  p  <  r-.  (p  and  r  are  Interchanged  when  p  .»  r  .) 
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Appendix  6 


I.  ( <S> (X)  -  1  f  4(m)alxc0I'P  if 

o 

where  for  m  •  1 

l(,)(x)  -  i  /(l-CM,).1’'"'*  d,  -  J0(x)  -  IJ,(X) 
for  m  ■  2 

l(2,(x)  -  I  /"(^cos^e'*008?  dtp  -  JQ(x)  +  I2J|  (x) 
and  for  m  >  2 

■  f  cos (m-l)cp  e,XC05%  -  Is-1  J-_,(x) 

o 

where  Jn(x)  is  the  Bessel  function  of  the  first  kind, 

A<">(y)  -  i  f"  e(n)e"Iyc0*6  Slnede 
n  o 

For  n  ■  I 

A^)(y)  -  i  cot  |  sine  e“,yc050  de  -  J0(y)  -  1 J,  (y) 
and  for  n  >  1 

A^(y)  ■  -  f  s In fn-1) ©  sine  e”,yc0*9  de 

TT 


III.  To  evaluate  the  Interrand  of  the  x-iniegral  of  Eq.  (28)  at  the 
singularity  X  ■  0  (see  Appendix  C)  ft  Is  necessary  to  define  the 
following  functions 
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Appendix  C 

Evaluation  of  the  integrand  of  the  X-integral  (Equation  28) 
at  the  Singularity  x  »  0 

The  integral  of  Eq.  (28)  is 


(C-l) 


where 


g(x)  “  l(n(IX+a£Nlp)  Km(lx*tajeNlr)  B--(x)e  3  for  p  :  r 
B--(x)  "  (flX+a2-*N+  («X+a2iN+  “Y)l^((q-  £)e£)  A^((q-  £)e£) 


m  •  q  + 

By  L'Hospital's  rule  the  Integrand  a.  \  ■  0  becomes 

1  im  g(x)  -  g(-x)  „  rag,(xl  _  ^g(-x)~| 

X-0  ,  l  bx 


(C-2) 


It  is  obvious  that 


[lm(IX+a(eNlp)Km(IX+a;/Nlr)]^  ^  *  [l  JI-X+a^IpjKjl-X+aiiNIr)]^  ^ 
e  X®0  *  e  a  ]'~n 


rirp^ 


R-1509 


FislAl .  aaiihl  1  o  2t  ls.  (ik)  I  b-  “(0) 

L  ax  ax  \mQ  a  '  'ml^Q  Dm,nv  ' 

+  (,K)mLH£ - — 1 


+  Bs,;<°> 


a[»m(IX+aANIp)Km(lX-faiNlr)] 


Sr>m(l-X+ajiNlp)Km(l-X+aANIr)]  | 


(C-3) 


Here  (lK)m|  -  lm ( I aXN I p)Km (1  ajiN I r)  for  p  £  r 

x*o 


(C-4) 


•s8S.n(-X) 


+8B5,nW 


-  ,(2,Vi+4)l("')(q<)A(n'>(^) 

x-o  r  P 


+  |(a2*N+  *Tjr)  (a2XN+  |m)  (qe£)A(n)  (q©g)+egl  (m)  (qeJ)A,(n)  (q©g)] 

?  c. 


and  |{">(x)  and  a|^(x)  are  as  defined  In  Appendix  B. 

The  third  term  of  (C-3)  Is  treated  as  follows: 
a)  For  X  ®  0+  and  ajlN  >  0 

!m(lX+ajiNlp)Km(lX+ajeNlr)  -  I m ( (X+a jeN) p) Km ( (x+a j&N) r) 


(C-5) 


lm(l-X+ajeNlp)Km(l-X+aXNIr)  =  I  r(a*N-x)p)K  ((a,eN-x)r) 


so  that  the  third  term  of  (C-3)  becomes 


“S.rW 


a[lm(g*l«ANI)p)Km((x+leWI)r)] 


■I  (for  p  £  r) 
'\*0 


C-2 
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'  28m.n-<°’{f  Km(IWNIr)L'm-l('aiNI',)  +  <la""l>>] 

-  y  lm(laXNIp)[Km-l  (laj&NIr)  +  K^,  (lajfcNI  r)]  (C-6) 

(Note  that  for  p  a  r,  p  and  r  are  Interchanged  In  Eqs.  C-3»  C-4  and  C-6.) 
b)  For  \  *»  0+  and  a^N  <  0 

lm(IX+aANIp)Km(lX+ajCNIr)  -  lm< ( I aXN I  -x) p)Km( ( I ajCN I -X) r) 


and 


lm(l-X+aiNlp)Km(l-X+ajWJIr)  -  lm((1aJeN!+x)p)Km((laJeNI+X)r) 


m'  - "  '  m’ 

The  third  term  of  C-3  then  becomes 


■  2B" -(0)  T*j}\  forpSr 

m*n  ax  *x-o 

Therefore  Eq.  (C-3)  con  be  written  as 


(C-7) 


rfr9,M. 

L  ax 


X"0 


21  (IK)  B-  -(0) 
3  m  a  in  f  n 


+  2(1 K) 


m 


x-o 


X-0 

M  I  a(IK)m 

±  2B-  -(0) - ® 

I  m,nN  ' 


m,n 

ax 


ax 


\*o 


where  ('K)I  Is  given  In  (C-'*) 


m 


X*0 


BSiS(0)  -  (a2iN+ 


as--  (x), 

f* " ""  Is  given  In  (C-5) 

x-o 


ax 

a(iK) 


m 


ax 


|  Is  given  In  (C-6) 
X**0 


and  the  upper  sign  Is  taken  when  i  >  0  and  the  lower  sign  when  l  <  0. 


C-3 


(C-8) 
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that 


When  !  »  m  ■  q  »  0,  by  the  limiting  process,  It  Is  easily  shown 


I im  B-  -(0)  -  llm  t -  0 
a  r\  m  f  n  a 

Sb  0  IbO 


a('K)  , 
llmB--(O)  o| 


,  ..  l 

l_0  m»n '  '  dX  V  i 


lim  ,  , 

!-0  urx'oi  n  -m,n 
X“0 


B-  -(0)  1  im  1  log!  -  0 

ir*  0 


llm  (IK) 


!-.  0  v"\e:0  dX 

Hence  when  !  a  m  ■  q  ™  0 


-  llm  (!+!  )  log! 

\=0  !r*0 


llm  ,g.(X). g,(,-x)  0 

X  0  x 

When  !  =  0  but  m  «  q  ?  0,  it  Is  easily  shown  that 


(C-9) 


1 


<*> 


Iml 


lim  (IK)  I 

!-*0  "Vo 


—L.  (l) 

2 Iml  V 


Iml 


I  Im 


adK) 


m 


!-*0  X»0 

Hence  for  !®0,  m  °  q  J*  0 


iim  ajxhaj^xi  2  jl 

X-0  X 


2  jiS  (IK>m 


for  o  ^  r 


for  p  i  r 


U 


*  j  I  ^  (qo£)A^  (qo£)[J  ^  -f-j  +  am  (“J  +  -y) 
1  r  D  r  o 


4t  l»J'  P  K)*®  KK< (S)  Kru{n'>  wfl 

ar  o 


(C- 10) 


C-4 


. . .  ‘irptFiff  in,,# 
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When  q  -  0  and  Eq.  (36b)  is  used  for  the  kernel  functions,  the  value 
of  the  Integrand  at  X  -  0  Is  zero  for  m  -  0.  For  m  ±  0  It  can  be  easily 
shown  that  the  integrand  at  X  ■  0  when  q  *  0,  m  ■  IN  Is 

1  T-(*2+  ^j)  (»2+  -y)lm(amp)Km(«mr) 

•  {ial  (m)  (0)A(n)  (0)  -  SJ  lfS)  (0)A(S)(0)  +  eg  I (S> {0)A,(S}  (0)}(C-l I) 


R-1509 


f 


I 


i 


1: 


t. 


The  Integral  I  Is 


r+26 


' '  -U  S*  dp '  ?  Lz6  C9>‘r)  +  69']d5 


r+2^ 


92  rr+2*  dP  %  /+26 

.2  “_7-  p-r  +  o’  J  .  . 2 

6  r-26  ^  r-26  (p-r) 


iL 


(P-4) 


where 


g2 


M,  +  M 

5 

H2  +  H4 

41 

• 

31 

• 

2(HS  - 

M,) 

(m4  - 

m2) 

41 

« 

31 

• 

-(M,  + 

+  4<M2  + 

m4) 

<♦1 

• 

31 

• 

2(M,  - 

M5) 

4(M2  " 

h4) 

41 

31 

The  Cauchy  principal  value  Is  to  be  taken  of  the  second  Integral: 
/+2fi  .. 

(  -i£- .  0 

r-26  P-r 

The  last  Integral  Is  the  Hadamard  type  whose  finite  contribution  Is 

r+26  A*  1  t 

r  .dP  _  Mm 


lim  | 
e-0  L 


r-€ 

.r+2fi  r+e~ 

r-26 

♦  /  +  J  J 

r+c  r-e J 

JL 

-1  r-e ,  r+26,  r+s 

P-r 

Jr-.?6,  r+e,  r-e 

I 

6 
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Therefore 


1 

.  1 

r 

-,r+26 

93 

1 

1 

"  T 

L69i 

*  P 

"T 

6 

V-26 

_l 

t 

r  ? 
L46  ' 

“i  9o 

siJ  -4' 
0 

'  ? (49' ' 

or 

i  < 

'?i 

■  5 

-  T2 

(H5-H|)  ■ 

■  >4 

■h2)] 

where 

Hi 

■  -863  * 

K'i,S)(p 

*  r  - 

26) 

h2 

"  -6 3  * 

K(5,K)(p 

*  r  - 

6) 

\ 

“  63  • 

k(5-S>(p 

“  r  + 

6) 

MS 

■  85*  * 

K(«.n-)(p 

*-  r  + 

26) 

and 


6  "  p/2 


(D-5) 


i 

1 

i 

1 

\ 


i 


D-3 


